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We study how the angular resolution of LISA for merging massive black-hole binaries would be 
improved if we observe multiple gravitational wave "images" due to strong gravitational lensing. 
The correlation between fitting parameters is reduced by the additional information of the second 
image which significantly reduces the error box on the sky. This improvement would be very helpful 
for identifying the host galaxy of a binary. The angular resolution expected with multiple detectors 
is also discussed. 

I. INTRODUCTION 



Black holes are a very fascinating prediction of general relativity. Recent observations indicate that massive black 
holes (MBHs) exist at the centers of many galaxies and their masses have a strong correlation with the properties 
of their host galaxies Binary MBHs would be formed by the merger of galaxies that harbor MBHs 2]. The 
coalescence of such a binary is the most energetic event in the Universe. It finally releases enormous energy (several 
orders smaller than c 5 /G ~ 4 x 10 59 erg/sec) over a characteristic time scale GM/c 3 ~ 5(M/10 6 Mo)sec (M: mass 
of the system) . The Laser Interferometer Space Antenna (LISA) can observe gravitational waves from MBH binaries 
with high signal to noise ratio (SNR) 0. In the case of equal mass binaries the expected SNR at a given distance is 
maximum for redshifted masses ~ 1O 6 M and is > 10 3 even at cosmological distances By fitting the gravitational 
waves from a MBH binary with templates prepared in advance we can make a stringent test of general relativity 
and obtain various important data about the black holes (e.g. mass, spin). Therefore gravitational waves from MBH 



binaries are, by themselves, a very interesting signal for relativistic astrophysics 0,0. 

If we can identify their host galaxies, it would also have a significant impact in astrophysics and cosmology. By 
detailed follow-up observations using electro-magnetic waves at various frequencies we might discover absorbing prop- 
erties of the galaxies or obtain fundamental clues which will elucidate the evolution of MBHs. In principle, the angular 
direction of coalescing MBH binaries could be estimated by analyzing data from LISA. Due to annual revolution and 
<*S rotation of LISA the signal from a MBH binary shows modulation that depends strongly on the ski position of the 
binary. For equal mass MBH binaries at z ~ 1 the angular resolution of LISA would be typically ~ 10 _4 sr for 
redshifted masses ~ 10 6 M Q 0. The resolution becomes worse for high redshift MBH binaries. The number of bright 
galaxies is ~ 10 4 in one square degree (~ 3 x 10~ 4 sr). Therefore the angular resolution of LISA is not sufficient 
to fully specify the host galaxy of a cosmological MBH binary 0|. Thus an improvement in the angular resolution 
for a MBH binary would be crucially helpful for selection of candidates for its host galaxy using electro-magnetic 
• i-H . telescopes. Possible mechanisms providing such an improvement are worth exploring. 

' In this context it has been discussed that higher gravitational wave harmonics of ||| , precession of the orbital plane 
in a rapidly spinning system 7], or effects caused by LISA's finite arm-length of @ could become important in some 
situations. In this paper we show that the angular resolution for a MBH binary could be significantly improved, if 
we observed multiple "images" of the gravitational waves due to strong gravitational lensing. Our point is not the 
amplification of the signal but observation of the same source at different epochs. The intrinsic binary parameters 
such as the chirp mass can be determined very accurately from gravitational waves, so it would not be a hard task to 
identify two lensed signals from the same merging MBH binary. 

If we assume that MBH binaries coalesce shortly after the merger of their host galaxies less than the Hubble time, 
the hierarchical model of structure formation suggests that the coalescence rate of MBH binaries could be as high as 
~ lOOyr , and its distribution might be much higher at high redshift than low redshift z < 1 |(g. The probability r 
of strong gravitational lensing also depends strongly on the redshift of the source. At z ~ 1 the probability is only 
t < 1CU 3 , but it could reach several percent at z ~ 5 [Toj . Thus we might obtain multiple images (more precisely the 
chirping gravitational waves) of distant MBH binaries. In this paper it is shown that the additional images caused 
by lensing with time delay AT > O.lyr would greatly decrease the correlations between fitting parameters of the 
templates, and the angular resolution of a MBH binary could be dramatically improved. 

This paper is organized as follows. In Sec II we briefly discuss signal analysis for a single image (II A) and for 
multiple images (IIB). Various numerical results are presented in Sec III. We also describe the angular resolution 
obtained with multiple detectors that are widely separated (IIIC). Our study is summarized in Sec IV. 
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II. ANALYSIS 
A. Observation of a chirping gravitational wave 

The three spacecrafts of LISA maintain a nearly equilateral triangle configuration and move annually around the 
Sun with the orientation of the triangle changing. As discussed below, the direction of a gravitational wave source 
can be estimated from the wave signal which is affected by (i) the Doppler effect due to the revolution of the detectors 
around the Sun and (ii) the amplitude modulation due to the rotation of the triangle and the increase of the wave 
frequency (chirp signal) 0, 3- For the response of LISA to gravitational waves we follow JTl| for the time-delay 
interferometry (TDI) analysis that cancels the laser frequency fluctuations in the data streams [lj. 

We label three spacecrafts as 1, 2, and 3 and use notation Li (i = 1, 2, 3) to denote the length of the arm opposite 
spacecraft i. Information about the binaries is extracted using three data streams A(t),E(t) and T(t) that form 
an orthogonal basis for the laser-noise canceling combinations. These three data streams are constructed from six 
(one-way) relative frequency fluctuations y%j(t) — 1,2,3, i ^ j) at different times (see |ll) for explicit expressions 
of A, E, T in terms of yij{t)). The signal is measured at spacecraft j and transmitted from spacecraft k(^= i, =/= j) 
along the arm L^. 

Our analysis is similar to (see also |l3l IT3, H^ L However, as in § we do not use the long wave approximation 
which is valid only for incoming gravitational waves larger than the arm-length of LISA L = 5.0 x 10 6 km, corresponding 
to / = c/2ttL ~ 0.01Hz. The study in Q is somewhat incomplete, as the TDI analysis was not included properly. 
This point is also remedied in this paper. For simplicity we assume that the arm-lengths Li are equal. 

First we briefly discuss the gravitational waveform from a MBH binary with a circular orbit. In this paper we only 
consider the in-spiral waveforms. For larger MBH binaries the contribution of the merger or the ring-down waveforms 
could be important for SNR E3- We use the stationary phase approximation and the restricted post-Newtonian 
approach with 1.5-PN phase |l6j | that is given in Fourier space as follows 

*(/) = 2nft c -<t> c -\ + \{^Gc-"M z f)-^[l + f (H + + (4/? - (2.1) 

where t c and <f> c are integration constants, and the former contains information about the coalescence time measured 
at the Sun not at the detector. The difference between the coalescence times is included in ea. H2.3f) . The parameter 
[3 is the spin-orbit coupling coefficient (we put its true value (3 = 0), and fi z , M z M z are the reduced mass, the total 
mass and the chirp mass respectively. All of the mass parameters are multiplied by the factor (1+z) (z: redshift of the 

binary) with the suffix z. The total mass M z is given by other two masses as M z = M.V 2 ' \i z 3 ^ 2 - The post-Newtonian 

expansion parameter x = 0(v 2 /c 2 ) is defined as x = {Gnc~ 3 M z f} 2 ^ 3 . 

For simplicity we first consider gravitational waves propagating in a homogeneous background. The effects of lensing 
will be discussed later. The incoming gravitational wave from a binary is decomposed into plus and cross polarization 
modes. It is convenient to use the principle axes (p, q) defined by the direction fl s and orientation S~2; of the binary as 
p = ft s x 0;/|O s x 1~2;| and q = £l s x p. The plus mode has the polarization tensor e^ b = p a Pb ~ <?a<7b and the cross 
mode has e* fc = p a % + qaPb- The relative frequency shift, for example, 2/31 (£), due to these two modes is expressed as 

2/31 (*) = \(vf) m (cos2Vi 2 A + + iS in2^ 12 A x ) (1 - cos# 12 ) [U{t, 1) - U(t - r, 2)] , (2.2) 

where r is the time r = L/c ~ 17sec and the function 

U(t, j) = exp[-2nif(t + Xj ■ H s /c)] (2.3) 

describes the phase of the incoming wave at the detector j (with position Xj measured from the Sun) at time t. In 
Ea. (|2.2|) we have neglected the very small relative motion of the detectors over the time scale r and also the time 
variation of the frequency over the time scale lAU/c = 500sec. The coefficients A + and A x are given in terms of 
an amplitude A as A + = A(\ + cos 2 /) and A x = 2/1 cos/ in which cos / = ft s ■ fl; is the cosine of the inclination 
angle. In Ea. i|2.2[l 9jk is the angle between the direction of the binary f2 s and the arm Xj — Xk, and the principle 
polarization angle ipjk is given by tan^fe = [(xj — Xk) ■ q]/[(xj — Xk) ■ p\- The angles 9jk and change with the 
motion of the spacecrafts. 

The Doppler effects are included in the factor U(t, j). The TDI data streams at time t are constructed from linear 
combinations of t/y at different times, e.g. t, t — r or t— 2r. Due to the finite arm-length, the response to gravitational 
waves shows a complicated behavior that depends strongly on the angular position of the MBH binary. The effect of 
finite arm-length on the estimation of the source direction could become important around ~ 10 5 + IQ^Mq in the case 
of equal mass binaries |8|. In the long wave limit /t< 1 the above expression l|2.2|l is simplified. The analysis by 
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Cutler |j| is essentially in this limit, and the amplitude modulation and phase modulation from Doppler effects can 
be analyzed separately. As a comparison we also examine the angular resolution of MBH binaries using this method 
to describe the response of LISA to gravitational waves. 

The intrinsic detector noise for a measurement of the fluctuations j/y has two main components that are relevant 
for the combinations A, E, T. They are the proof mass noise and the optical path noise. We use the spectra 

gproof^mass^ = 3 5 x 1Q -48 (-^-) Hz -1 , (2.4) 
gopUcal-path^ = l8xW -37 fj_\ j^-l ^ 

It is a simple task to obtain the noise for the three combinations A, E, T. But we note that the detector's noise below 
O.lmHz is somewhat controversial (see e.g. [l4l[r^ |). This frequency region strongly affects our results for larger MBH 
binaries. 

We also include the Galactic binary confusion noise [l8| using the following simple approximation 

^ C (/) = | 4X10_37 (^)^ HZ_1 (^ 1(r2 ' 7Hz ) (2.6) 
I (/ > l(r 2 - 7 Hz) 

given for one year observation. This confusion noise is about 10 times larger (for h) than the detector noise at lmHz 
(see figure 5 of 01 )■ For simplicity we neglect the anisotropies of the Galactic confusion noise in this paper. 

We integrate each chirping gravitational wave for lyr before coalescence up to the cut-off frequency f ma x when the 
binary separation becomes r = 6GM z /c 2 (l + z), corresponding to a frequency 

fma X =0M -^-) Hz. (2.7) 



Vl0 5 M Q 

In contrast the wave frequency / observed at Tyr before the coalescence time t c is given approximately by 

M \ - 5 /« / T N - 3/8 



/ = 1.9 x 1(T 4 \ %r Hz. (2.8) 

The signal to noise ratio SNR of a detection is determined by the amplitudes of the modes A, E and T and by the 
noise spectra SU(/), SE{f) and Sx(f) which are constructed from the detector noises l|2.4|l l|2.5|l and the confusion 
noise (|2.6|) . We calculate the signal to noise ratio of the in-spiral waves using 0, |2(J 

snr 2 = 4 v t"W!l df , M 



B=A,E,T 



fi 



where f\ is the initial frequency for signal integration. In the present study the number of fitting parameters {a{\ 
for each chirping gravitational wave is 10. They are {a-i} — {M. z ,H z ,l3,t c ,(j) c ,A,6 s ,<p s ,6i,<pi}. Angular variables 
($s; <Ps) represent the direction of the binary f2 s in a fixed frame around the Sun, and (6*/, 4>i) represent its orientation 
f2/. The variance in the parameter estimation errors Acti are evaluated using the Fisher information matrix as 
(AetiAaj) — T^ 1 where T^- is given as [2(J 

r„-4 E t" (2.0, 



B=A,E,T 



Following Cutler we use the notation AQ. S = 2i:smd s y (A8^) {AipD — {A9 s A<f> s ) 2 for the angular resolution of a 
binary on the sky. 

Our results for the angular resolution AO s change only slightly for masses > 10 4 M Q when we switch the integration 
time from lyr to 10 yr. We put the MBH binaries at z = 3 in a universe with cosmological parameters f^o = 0.3, 
Ao = 0.7 (spatially flat universe) and Hq = 71km/sec/Mpc. In this cosmological background the luminosity distance 
d,L becomes dh = 25Gpc at z = 3 and — 6.5Gpc at z — 1. 
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B. Analysis for multiple images by lensing 

We discuss parameter fitting of multiple images created by strong gravitational lenses. As we show below the 
improvement in the angular resolution Af2 s due to lensed multiple images could be very effective for time delays 
AT > O.lyr. On the other hand a time delay larger than the operation period of LISA (nominally ~ 3 yrs, but possibly 
~10 yrs) is irrelevant for our analysis. Thus our main targets are multiple images with time delay O.lyr < AT < lOyr. 
LISA has good sensitivity to gravitational waves with frequency 10 _4 Hz < / < 10 _1 Hz which is much higher than 
the inverse of the relevant time delay (AT) -1 < 3 x 10 _7 Hz. We note also that the latter frequency is smaller than 
/ given by ea. l|2.8|l for T = lyr and M. z < 10 s M Q . Therefore the geometric optics approximation is valid in our 
analysis and the structure of phase is not changed by the lensing 0,H2|. We take the same parameters Ai z , fj, z , 4> c , ft 
for the phases ea. l|2.1ll of all the multiple images from a source, but the coalescence time t c is, of course, different. In 
some cases the substructure of the lensing galaxy might complicate the problem |2^| . 

In the case of a homogeneous and isotropic universe the wave amplitude A is given simply in terms of the luminosity 
distance g?l and the chirp mass M z by (2J] 

Mi) 

The above relation can be used to determine cIl from the observed gravitational wave and also to specify the redshift 
of the binary z by inversion of the (Ll-z relation if the cosmological parameters are known accurately |14L |24|. Note 
that in general the chirp mass M z is determined more accurately than the amplitude A using the time evolution of 
the gravitational wave phase. 

For signals affected by strong lensing we fit different values Ai for the amplitude of each image. We cannot apply 
the above inversion to estimate z, but this is not a serious problem for identifying high-z MBH binaries, considering 
the accumulated effects of weak lensing. At z ~ 3 the rms fluctuations of the amplification by weak lensing could 
become ^0.2 and hamper the application of the above inversion even for a single image [l5ll25j . 

Here we discuss fitting of the angular variables (9 S , 4> s ,9i,<j)i). First we evaluate the typical image separation for a 
lensed source. As a model we use the singular isothermal sphere approximation for the density profile of the lensing 
galaxy. In this model we have the following explicit relation between the image separation angle 9 and the time delay 
AT |2U 

AT = + ZL)y = o ( 4V (**\ (^-) (±±*l) (JL) vr , (2 . 12) 

2c D LS y L,J \5"J \DlsJ \lAGpcJ \ 2 J \0.5J y ' V ' 

where D = + z) 2 is the angular diameter distance {e.g. Dol- between the observer and the lens), zl is the 

redshift of the lens, and y is the normalized impact parameter with y < 1 corresponding to the occurrence of multiple 
images in the isothermal sphere approximation. The image separation 5 ) for the relevant time delay AT ~ lyr is 
much smaller than the angular resolution of the source direction for cosmological MBH binaries, which will be shown 
later. This means that we can use the same fitting parameters (9 S , 4> s ) for the direction of the binary for each image. 

The polarization tensor e a & of a gravitational wave is parallel transported along a null geodesic [27J. It is easy to 
confirm that the directions of the tetrad vectors on a lensed null geodesies change only on the order of the image 
separation, which is generally much smaller than the resolution of the orientation of a cosmological MBH binary 
(9i,(fii). This means that the polarization properties of the multiple images cannot be distinguished observationally 
and the same fitting parameters {9i, 4>i) can also be used for the orientation of each image. 

Let us briefly summarize the number of fitting parameters. For n- multiple images from the same source we can take 
eight common parameters {M z , fJ, z , f3, <j> c , 9 S , (/)„, 9i, 0;}, but fit two different ones {Ai, t C i\ for each image (i = 1, n). 
Thus the total number of fitting parameters becomes 8 + 2n. In the following analysis we fix the true value of each 
wave amplitude A (both un- lensed and lensed) by A = Aq. This would give conservative results for the magnitude 
of the parameter estimation errors in the case of lensed signals. With our prescription the SNR depends on the total 
number of images n as SNR oc n 1 ^ 2 . 

III. ANGULAR RESOLUTIONS FOR MBH BINARIES 
A. Basic analysis 



We calculate the parameter estimation errors using (i) the TDI method described in Sec II and (ii) the method of 
Cutler |4| which uses the long-wave approximation. The effective noise curves are somewhat different for these two 
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approaches. In this section various averaged quantities, such as the SNR and the angular resolution Af2, are evaluated 
by taking geometrical averages for 100 MBH binaries at z = 3 with random directions (9 S , <p s ) and orientations (#;, <j>{). 
For simplicity we only study equal mass binaries. The time delay is fixed at AT = l/3yr unless otherwise stated. 

In Figure 1 the averaged the SNR is presented for redshifted masses 4 x 10 3 — 4 x 10 8 M Q (true masses 10 3 — 10 8 M Q ). 
Hereafter we mainly use the redshifted mass m z = m\ z = mi z to show the mass dependence. When we decrease the 
mass from m z ~ 1O 8 M0, the SNR becomes a maximum around m z ~ 1O 6 M0. This is because the binary confusion 
noise disappears around / ~ 10~ 27 Hz in our models and this frequency corresponds to the coalescence frequency 
fmax of a MBH binary with mass m z ~ 10 6 M© as given by ea.l|2"77Jl. The SNR for two lensed images simply increases 
by a factor of ~ 1.5 ~ y/2 as expected. 

In Figure 2 we show the averaged angular resolution as a function of redshifted mass m z . For single images 
the averaged resolution Af^i shows weak dependence on the mass for the TDI method, and is nearly constant at 
~2x 10~ 3 sr for MBH binaries with masses between m z ~ 10 5 M Q and 10 8 M Q at z = 3. This is a remarkable contrast 
to Figure 1 for SNR which shows a steep rise around m z ~ 10 5 — 1O 6 M . Using the simple method of Cutler we find 
that the sky positions of a MBH binaries with m z > 10 5 M Q is mainly estimated from the amplitude modulation due 
to rotation of the detectors. The Doppler phase modulation has a contribution for m z < 1O 5 M . 

When a second image is added by the time delay of lensing, the situation changes drastically. The angular resolution 
Af^2 obtained from two images with m z ~4x 10 5 Mq is improved by more than two orders of magnitude, compared 
to a single image Af2i. The ratio Aili/Af^ decreases for both smaller and larger masses. It becomes ~ 10 at 
m z ~ 1O 4 M0 and ~ 5 at m z ~ 1O 8 M . This mass dependence is discussed later. We also found that the parameter 
estimation errors for the intrinsic binary parameters such as the chirp mass would be changed by only a factor of 2 
or so by the second image. 

In Figure 3 we show histograms of the angular resolutions Afii and Af^ in our 100 realizations of MBH binaries 
with redshifted masses m z = 4 x 1O 3 M ,4 x 10 5 Mq and 4 x 1O 7 M . The impact of lensing is also apparent in these 
figures. 

So far we have fixed the time delay at AT = l/3yr. In Figure 4 we present the ratio Atti/A^ as a function of the 
time delay AT. Due to the periodicity in the configuration of LISA, the results obtained for ATyr are the same as 
those for AT + Nyr with ./V an integer. As our results are given in the form of a ratio AQi/ Af^, they do not depend 
on the distance or redshift to the source as long as we use the redshifted masses. The factor Af2i/Afi 2 depends only 
weakly on the time delay for AT > O.lyr. 

B. Third image 

We have also calculated the angular resolution Af^ in the case that we could observe a total of three images. We 
set the two time delays relative to the first image as AT = 1/3 yr and AT = 2/3 yr, and fixed the amplitudes of the 
three images by A Q given in ea. l(CTTl . The ratios ASI1/AO3 become ~ 40 (m z = 4 x 10 3 M©), ~ 800 (4 x 1O 5 M ) 
and ~ 22 (4 x 1O 7 M0). For the first two images analyzed in the previous subsection the ratios ASli/Af22 are ~ 13 
(4 x 1O 3 M ), ~ 300 (4 x 1O 5 M ) and - 9 (4 x 1O 7 M ). Therefore the effect of a third image is not so drastic as that 
of the second one. 



C. Multiple detectors 

We have shown that the observation of multiple images due to lensing could significantly improve the angular 
resolution A£l. But we cannot bet on this passive effect for most merging MBH binaries, considering the probability. 
In an usual situation we will only detect a single image. Using multiple detectors is one positive method for improving 
the angular resolution. In this case we can take the same fitting parameters also for the amplitude A and the 
coalescence time t c (at the Sun). The former is important as the amplitude A and the angular variables will correlate 
strongly (see ea. H2.20 ). The latter means that we can take advantage of the time delay between the two detectors 
which is closely related to the direction to the binary. 

As a simple extension of our analysis for lensing, we calculate the angular resolution AQu with two detectors 
that have the same specification as LISA. The two detectors are on the Earth orbit around the Sun as for LISA, 
but we set the angle between them at 2-7r/3 corresponding to l/3yr of orbital time. Thus their distance is fixed at 
lAUx 2 x sin(7r/3) =1.73AU, and their orientations are always different from each other. For simplicity we assume that 
the noises in the data streams measured by the two detectors are independent. This would be a valid approximation 
for the detector noise, but might not be for the binary confusion noise. We evaluate the Fisher matrix using an 
extension of eq.(2.10) with six data streams; (A,E,T) modes both for the two detectors. 
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The ratio of angular resolutions for a single LISA Afij and for two LISAs Attn becomes Afii/ Attn ~ 10 2 
(4 x 10 3 M Q ), ~ 6 x 10 3 (4 x 10 5 M Q ) and ~ 30 (4 x 1O 7 M ). Therefore even for binaries at z = 3 the angular 
position of a merging MBH binary can be determined with error box ~ 10~ 6 sr, namely (several arc- minutes) 2 . This 
area on the celestial sphere is close to the resolution obtained from Gamma-Ray-Bursts for determining their host 
galaxies. When another detector (a total of three) is added with positions characterized by l/3yr and 2/3yr, the ratio 
AQi/Afijxr becomes ~ 400, ~4x 10 4 and ~ 60 respectively. 

D. Effective time 

As shown above, the angular resolution Att could be improved significantly for m z = 10 5 ~ 1O 6 M0. Here we analyze 
this mass dependence. Roughly speaking, the SNR characterizes the magnitude of the Fisher matrix Tij as we can 
understand from eas. (|2.9|) and (|2.10l) . When the correlation (degeneracy) between the fitting parameters is large, the 
parameter estimation errors such as Att also become large. This correlation depends strongly on the configuration of 
LISA around the epoch when the SNR accumulates. In this situation some independent information e.g. added by 
the second image, has an important role in reducing the correlation and improving the angular resolution Att. This 
fact can be anticipated by comparing Figures 1 and 2. 

The angular directions of MBH binaries are mainly estimated from the amplitude modulation and the Doppler phase 
modulation caused by the revolution and rotation of LISA. A long duration observation is crucial for determining 
them well. Therefore the correlation between parameters including the direction would be related to the effective 
observational period, and a binary with a shorter observational period would have a larger ratio Atti/ Att,2- Let us 
define a time t e ff by weighting the time before coalescence t c — t by the amplitude of the wave h(f) and the noise 
curve Sh(f) as follows 

Here we formally wrote down the above expression (see ea. H2.90 ). The denominator is nothing but the square of the 
SNR. As the weight factor \h(f)\/ Sh(f) 1 ^ 2 does not have a very spiky structure, we can regard t e ff as an effective 
observational period of the binaries. For monochromatic sources we have t e ff = l/VSyi. In figure 5 we show the 
time t e fj (solid curve) and the ratio of the angular resolution Att 2 / AQj (dashed curve). The overall shapes of the 
curves are similar as expected. 

The mass dependence of t e ff can be understood as follows. For larger mass binaries with m z > 1O 5 M0, the SNR 
(denominator of ea. l|3.1|l ) comes mainly from the wave emitted close to the final coalescence. But these waves have 
a small contribution to the numerator. A binary with mass m z < 1O 6 M has coalescence frequency larger than 
10 - Hz, which is critical for the Galactic binary confusion noise. Therefore the effective observational time t e ft 
decrease significantly for m z < 10 6 Mq. The smaller mass (m z < 10 4 M Q ) binaries stay in the sensitive frequency 
region for a longer time, and the effective time t e ff increases. 

IV. SUMMARY 

In this paper we have studied how the angular resolution of LISA would be improved if we observe multiple images 
due to strong gravitational lensing. It is found that the error box on the sky could be typically 100 times smaller for 
redshifted mass m z = 10 5 ~ 1Q 6 M . The improvement in the angular resolution depends strongly on the masses of 
the MBH binaries. This mass dependence can be roughly explained by using an effective observational time for the 
binary. 

We have mainly discussed the effects of lensing on the matched filtering analysis of gravitational waves, but lensing 
would also be of great advantage in searching for the host galaxy of a MBH binary using electro-magnetic waves. 
Besides the fact that the target itself is lensed, additional information obtained from the gravitational waves, such as 
the time delay or the ratio of the amplification factors of the images, could be used to specify the host galaxy. 

We have also calculated the angular resolution expected for multiple detectors and found that the resolution could 
be improved by a factor of ~ 6 x 10 3 (m z ~ 4x 10 5 M Q ) compared to that from a single detector. If LISA could detect 
gravitational waves from MBH binaries, very exciting results might be obtained by launching another detector. 
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FIG. 1: Averaged SNR for equal mass binaries at z = 3. The thick curves are for the TDI method and the thin ones are for the 
method of Cutler. The results for the single (double) image are given by the dashed (solid) curves respectively. The horizontal 
axis represents the redshifted mass m z . 
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FIG. 2: Averaged angular resolution. The identification of each curve is the same as in figure 1. 
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FIG. 3: Histograms of the angular resolution Afii (open bars) and AQ 2 (gray bars). The mass m z is fixed at (a) m z = 4 x 10 3 , 
(b) 4 x 10 5 and (c) 4 x 10 7 . The results for Af2i are displaced by 1CP 0,1 . The total number of binaries is 100 for each mass. 
These results were obtained using the TDI method. 
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FIG. 4: Averaged ratio Afli/Afi2 as a function of the time delay AT. The three curves represent different masses. These 
results do not depend on the distance to the binaries. The TDI method was used again. 
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FIG. 5: The effective time t e ff (solid curve: left axis) and the ratio (Af2i/Afi2) (dotted curve: right axis) as functions of 
the redshifted mass m z . The TDI method was used here again. 



